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Abstract 

We report a preliminary measurement of the first and second moments of the hadronic mass dis- 
tributions in B — ► X c lv decays. The measurements are based on T(4S) — > BB events where the 
hadronic decay of one of the B mesons is fully reconstructed and a charged lepton from the decay 
of the other B meson is identified. The moments are presented for threshold lepton momenta 
ranging from 0.9 to 1.6 GeV. From the (M^) moments we determine the non-perturbative Heavy 
Quark Expansion (HQE) parameters, A and Ai. We combine the measured moments (M^) with 
earlier BABAR measurements of the semileptonic branching ratios and B lifetimes and perform a 
simultaneous fit to the HQE for the moments obtained for different threshold lepton momenta and 
the semileptonic decay width. This fit results in an improved value for the CKM matrix element 
\V cb \. 
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1 Introduction 



The heavy quark limit in QCD has become a very useful tool for relating inclusive -B-decay prop- 
erties, like the semileptonic branching fraction and moments of mass and lepton momentum distri- 
butions, to the charged current couplings, \V c b\ and \V u b\ ■ 

In this paper we report a measurement of the moments of the invariant mass distributions, 
(Mx) and (M x ), of the hadronic system recoiling against the charged lepton and the neutrino in 
semileptonic B decays. The moments are presented for threshold momenta of the charged leptons, 
ranging from 0.9 GeV to 1.6 GeV. This momentum range retains the sensitivity to the spectrum of 
hadronic final states, including higher mass resonant charm mesons and non-resonant D^ir states. 

A preliminary measurement of the threshold momentum dependence of (M x ) was presented 
last summer [1]. This earlier analysis followed a procedure first introduced by the CLEO Collabora- 
tion [2] in which the moments were derived from the branching fractions and average hadron mass 
distribution of various charm states. For this purpose, the mass distributions were decomposed 
into contributions from various charm meson resonant and non-resonant states. A limitation of 
this approach was that it required as input the knowledge of the mass distribution for the indi- 
vidual charm states, and this information as well as the branching ratios remain uncertain for the 
higher mass states. 

The new analysis presented here is much less dependent on this input. The moments are 
extracted directly from the measured Mx and M x distributions, taking into account corrections 
for the mass scale, the detection efficiency, and small residual backgrounds. 

Theoretical calculations of the second moment (M x ) have been carried out [3] using an Operator 
Product Expansion (OPE) in powers of the strong coupling constant a s {m\,) and in inverse powers 
of the B meson mass up to order a 2 s (3o and \jm\. (Here /?o = (33 — 2n/)/3 = 25/3 is the leading 
order QCD (3 function, and rif is the number of quark flavors accessible at this energy, i.e. nj = 4.) 
These expansions contain the non-perturbative parameters A (0(l/mg)), Ai, and A2 (0(l/m 2 B )). 
The moment measurements presented here extend over a large region of phase space and thus we 
can use them to extract A and Ai and to test the validity of the assumptions made in the OPE. 

By combining measurements of (M x ) for seven values of the momentum threshold with earlier 
BABAR measurements of the semileptonic branching ratios and B lifetimes we can further constrain 
the parameter Ai and the b quark mass, ml , in the 15 scheme [4] (where parameters like the quark 
mass are defined by reference to the mass of the T(15) state) and derive an improved measurement 
of the CKM matrix element \V c b\. Furthermore, we use the published data on moments of hadron 
mass-squared, lepton energy and photon energy distribution in semileptonic and rare radiative B 
decays, to test the consistency of the underlying theoretical framework. 

2 The BABAR Detector and Data Set 

The measurements presented here are based on a sample of about 89 million BB pairs collected at 
the T(4S) resonance by the BABAR detector [5] at the PEP-II asymmetric-energy e + e~ storage ring 
operating at SLAC. The detector consists of a five-layer silicon vertex tracker (SVT), a 40-layer drift 
chamber (DCH), a detector of internally reflected Cherenkov light (DIRC), and an electromagnetic 
calorimeter (EMC) assembled from 6580 CsI(Tl) crystals, all embedded in a solenoidal magnetic 
field of 1.5 T and surrounded by an instrumented flux return (IFR). 
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3 Analysis Method 



This measurement of the hadronic mass distribution in semileptonic B decays exploits the very 
large sample of B mesons recorded by the BABAR detector. The analysis relies on BB events in 
which one B meson decays hadronically and is fully reconstructed (B reco ) and the semileptonic 
decay of the other B meson (B reco n^ is identified by the presence of an electron or muon. While 
this approach results in a low overall event selection efficiency, it allows for the determination of the 
momentum, charge, and flavor of the B reco meson. Furthermore, the remaining tracks and photons 
in the event can be used to reconstruct the hadronic mass Mx of the semileptonic decay. The 
determination of Mx is substantially improved by application of a kinematic fit to the full event. 

3.1 Selection of Fully Reconstructed Hadronic B Decays, B — > DY 

To reconstruct a large sample of B mesons, hadronic decays of the type B reco — ► D^Y^ are 
selected. Here, D^*> refers to a charm meson and consists of hadrons with a total charge of ±1, 
composed of niir^ n^K® n^ir with ri\ + n-i < 5, 77,3 < 2, and 714 < 2. We reconstruct D*~ — > 

D°tt~ and D*° -► D°ir ,D°~f and the decays D~ -> K+tt-tt', K+tt-tt-tt , K s -k~ , i^Ti-vr , 
K°7r-7r-vr + ; and D° -» K + n~, K~^~ ir~ ir® , K + 7r-7r-vr + , K°ir + ir-. 

The kinematic consistency of the B reco candidates is checked with two variables, the beam 

energy-substituted mass ttt-es = v/s/4 — and the energy difference AE = Eb — \fs/2. Here 
y/s is the total energy in the T(4JS) center-of-mass frame, and ps and Eb denote the momentum 
and energy of the B reco candidate in the same frame. We require AE to be within three standard 
deviations from zero for each decay mode. 

For a given B reco decay mode, the purity is estimated as the fraction of signal events with tties 
> 5.27 GeV. For this analysis we only use decay modes for which the purity exceeds 40%. In events 
with more than one reconstructed B reco decay, we select the decay mode with the highest purity. 
For the case of two candidates being reconstructed in the same mode, the one with AE closest to 
zero is selected. 

3.2 Selection of Semileptonic Decays, B — ► Xiw 

Semileptonic B decays are identified by the presence of one and only one electron or muon above 
a minimum momentum p* min measured in the rest frame of the B reco n meson recoiling against the 
B reco . The lepton momentum threshold p* min is varied in the range of 0.9 — 1.6 GeV. The higher 
the lepton momentum, the lower are the residual backgrounds from secondary decays of charm 
particles, leptons, and from misidentified hadrons. Measurements with lower lepton momenta 
are more sensitive to the production of higher mass charm mesons, in the form of both resonant 
and non-resonant states. 

To improve the background rejection, we impose a condition on Qe, the charge of the lepton 
from B recoi i, and Qb, the charge of the 6-quark of B reco , namely QbQi < 0. This condition is 
fulfilled for primary leptons, except for B°B° events in which flavor mixing has occured. 

Electrons are identified using a likelihood-based algorithm, combining the track momentum 
with the energy, position, and shape of the shower measured in the EMC, the Cherenkov angle 
and the number of photons measured in the DIRC, and the specific energy loss in the DCH. The 
efficiency of the electron selection has been measured with radiative Bhabha events and has been 
corrected for the higher multiplicity of BB events using Monte Carlo (MC) simulations. 



9 



Muons are identified using a cut-based selection for minimum ionizing tracks, relying on informa- 
tion from the finely segmented instrumented flux return of the magnet. The number of interaction 
lengths traversed by the track, the spatial width of the observed signals, and the match between 
the IFR hits and the extrapolated charged track are used in the selection. The muon identification 
efficiency has been measured with events and two-photon production of /U + /^~ pairs. 

The misidentification probabilities for pions, kaons, and protons have been extracted from 
selected data control samples. They vary between 0.05% and 0.1% (1.0% and 3%) for the electron 
(muon) selection. 

To further reduce backgrounds we require the total charge of the event to be \Qtot\ — \QBreco + 
QBrecoiil < 1- We explicitly allow for a charge imbalance to reduce the dependence on the exact 
modeling of charged particle tracking in the Monte Carlo simulation, especially at low momenta, 
and the production of tracks from photon conversions. 

3.3 Reconstruction of the Hadronic Mass Mx 

The hadronic mass Mx in the decay B — > Xlv is reconstructed from charged tracks (ir^ and K^) 
and photons that are not associated with the B reco candidate or identified as leptons. Specifically, 
we select charged tracks with a minimum transverse momentum pt > 0.1 GeV and at least 12 
DCH hits. For photons we require an energy sum greater than 0.05 GeV deposited in three or 
more crystals. To suppress background from hadrons interacting in the calorimeter we require a 
minimum angular separation between the center of the shower and the nearest track impact point 
of 200mrad and a shower shape consistent with a photon (lateral moment LAT < 0.5, defined in 
reference [5] ). 

The measured four-momentum Px of the hadron X can be written as 

N ch 

Px = Y, p f h + Y, p ] (!) 

i=X j=l 

where P are four-momenta and the superscripts ch and 7 refer to the selected charged tracks and 
photons, respectively. Depending on particle identification the charged tracks are assigned either 
the or mass. 

The neutrino four- momentum P v is estimated from the missing momentum four- vector, P m i ss = 
-fr(4S) — -fBr-eco — Px — Pe, where all momenta are measured in the laboratory frame. The mea- 
sured missing momentum four- vector is an important indicator of the quality of the reconstruction 
of the total recoil system. Any particle that is undetected or poorly measured, any sizable en- 
ergy deposition due to charged hadron interactions or beam-generated background will impact the 
measurement of Mx and M^ niss = P^iss- To reduce the impact of these experimental effects, we 
impose the following criteria: E m i ss > 0.5 GeV, \p m iss\ > 0.5 GeV, and \E m i ss — \p m i ss \\ < 0.5 GeV. 
In Figure 1 the E m i ss , \p m iss\, and E m i ss — \p m i ss \ distributions are shown and compared to MC sim- 
ulations, normalized to the total number of events. The agreement reflects our understanding of the 
detector performance and improvements in the suppression of track and photon background which 
have been studied extensively. It should be noted that the branching fractions for the individual 
decay modes used in the MC can also contribute to differences between data and MC. However, the 
extraction method for the moment measurements is not sensitive to branching fraction variations. 

To further improve the resolution on the measurement of the hadronic mass Mx, we exploit 
the available kinematic information from the full event, namely the B reco and the B reco n candidate, 
by performing a kinematic fit with two constraints that imposes four-momentum conservation, the 
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Figure 1: Data - MC comparison of missing energy E m i ss (left), missing momentum \p m iss\ (middle) 
and E m i ss — \p m iss\ (fight) for p* min = 0.9 GeV. The MC predictions for the total distribution as 
well as the residual background are indicated. 



equality of the masses of the two B mesons, M recoi i = M reco , and forces M^ iss = 0. The fit takes 
into account event-by-event measurement errors for all individual particles and the four-vector of 
the measured missing momentum. The resulting Mx resolution is 350 MeV. 

3.4 Subtraction of B reco Background 

Fi gure 2 shows the m-ES distribution for the B reco candidates for all events for which the B reco n 
meson meets all selection criteria for a semileptonic decay. To estimate the combinatorial back- 
ground in the signal region (vies > 5.27 GeV), the measured uies distribution is fitted to a sum of 
a background function (introduced by the ARGUS Collaboration [6]) and a signal function (first 
used by the Crystal Ball Collaboration [7]). For toes > 5.27 GeV, we find a total of 7114 signal 
events above a background of 2102 events. In the analysis, the B reco background subtraction is 
carried out for several regions in Mx to account for changes in the signal-to-background ratio as a 
function of Mx- 

4 Determination of the Hadronic Mass Moments 

The distributions in Mx and M x , corrected for combinatorial B reco background, are shown in 
Figure 3 for two different lepton threshold momenta. The dominant contributions are from the 
lowest mass mesons, (D + , D ) and (D , D* + ), but there are clear indications of the higher mass 
charm meson states. The residual background, estimated from MC simulation, is primarily due to 
hadron misidentification and non-prompt leptons from semileptonic decays of and D s mesons, 
either from B°B° mixed events or produced in W~ — > cs fragmentation. 

To extract the first moments from these distributions, i.e., the unbiased mean values {Mx) 
and (M x ), we need to correct for effects that can distort these distributions and thereby affect 
the moments. Specifically, we need to calibrate the measurement of Mx and M x by relating the 
observed values to the generated values for MC simulated distributions. Furthermore, we need 
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Figure 2: XTie m^g distribution for B reco candidates in events for which the B reco u meets all 
selection criteria for a semileptonic decay. The lepton momentum threshold is p* m i n = 0.9 GeV. 



to reliably estimate and subtract the contribution to the moments from residual backgrounds and 
then correct the result for the detection efficiency. Thus we can express the true hadronic mass 
moment (M x ) as 

10 

I ]\/T n \ DATA f I Ti/rn \ MCbg 

(M x ) = calib ' calib x — (2) 



1 - f t io 

Jbg ^ Rf et (M% 



MCdet 
X,i/ calib 

where the first term describes the subtraction of the residual background and the second term spec- 
ifies the efficiency correction. Here (M x ) *"^f' e (with sample = DAT A, MCbg, MCtrue, MCdet) 
refer to the mass moments (calibrated and corrected for B reco background) extracted from samples 
of data, of MC-simulated background, and of ten MC-generated and -detected exclusive semilep- 
tonic decays, B — > X' l c iv. These decay modes, identified by the index i, include D, D*, four resonant 
D** and four non-resonant D^ir charm meson final states, summed over charged and neutral B 
mesons. The coefficient Ri is the relative fraction of a decay mode predicted by MC before detection 
efficiencies have been taken into account; Rf et are the corresponding fractions for detected decay 
modes. We denote by fb g the size of the residual background relative to the data. 

To study the calibration, efficiencies, and backgrounds we rely on Monte Carlo simulations. 
Decays to D*lv are modeled by HQET-derived form factors [8], decays to Div and D**lv are 
simulated using the ISGW2 model [9], while the non-resonant decays D^rriu are modeled according 
to the prescription by Goity and Roberts [10]. The relative branching ratios are adjusted to be 
consistent with current measurements. 

The calibration of the mass scale is performed separately for (Mx) and (M x ). The results of 
the procedure are illustrated in Figure 4. The average measured versus the average true value for 
bins in Mx and M x are determined from a MC simulation of decays B — > X c lv. The relation 
between the measured and true quantities can be approximated by a linear function in both cases. 
Though the calibrations are obtained from averages in bins of Mx and M x , they can be applied 
on an event-by-event basis using the slope and offset of the linear functions obtained from fits. 
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Figure 3: Mx (left) and M x (right) distributions after subtraction of the B reco background, with 
Pmin = 0-9 GeV (top) and Pm in = 1.5 GeV (bottom). The dominant contributions from D and D* 
mesons (histogram) and the residual background (shaded area) as estimated from MC are indicated. 

To verify this, the calibrated M x are averaged for each of the ten decay modes. As can be 
seen in the upper halves of Figures 4a and b, the calibrated moments (M x ) ^^ et reproduce the 
underlying true moments (M x ^ MCtrue for each decay mode. This result indicates that there is no 
significant mass dependent or decay-type dependent bias that is unaccounted for. The calibration 
is largely independent of the decay model, and the impact of decays which are not included in the 
MC simulation (or which are mediated by other processes) can be treated adequately. 
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Figure 4: Results of the calibration procedure for events with a lepton threshold of p min = 0.9 GeV , 
a) for (Mx) and b) for (M x ). The calibration data and fit results are shown in the lower half 
(round data points and dashed line), the verification in the upper half (triangles and solid line) of 
the figures. 

Detailed studies show that the slope and offset of the calibration curves vary slightly as a 
function of the multiplicity of the hadron system and E m i ss — \p m iss\- Thus instead of one universal 
calibration curve for all data, we split the data in three bins in multiplicity and three bins in 
Emiss — \i>miss\, and derive a total of nine calibration curves, one for each subsample. 

Decays to higher mass final states usually correspond to higher multiplicities in the decay. 
Therefore, the event selection criteria are expected to have a stronger impact on high mass states. 
In addition, the different decay modes have different spin configurations and thus different angular 
distributions. The second term in Equation 2 accounts for these effects. It is determined by MC 
simulation of the individual semileptonic decays B — > X % c tv. 

The hadronic mass moments (Mx) and (M x ) obtained after background subtraction and mass 
calibration are presented in Figure 5 as a function of p min with statistical and systematic errors. 
The results are also summarized in Tables 1 and 2. As expected, the values for the hadronic 
mass moments increase for lower p min ■ Similarly, the systematic errors associated with background 
subtraction increase with lower p min , while most of the detector related systematic errors are 
independent oi p min . This results in a systematic error which is of comparable size to the statistical 
error for most of the p^ in range. It should be noted that the statistical error of the measurements 
is almost constant as a function of p min because the width of the M x distribution decreases as 
p m i n increases. The correlation coefficients for the different p m i n measurements of {Mx) vary from 
~ 90% between neighboring points to ~ 55% for the two extreme p min values. These correlation 
coefficients are slightly smaller for (M x ). 
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Figure 5: Measured moments (Mx) (left) and {M x ) (right) for different lepton threshold momenta, 
p* min ■ The bars indicate the statistical and the total error. The errors of the individual measurements 
are highly correlated. 



Table 1: Preliminary results for (Mix) for different values of p m i n , with statistical and total 
systematic errors. The last four columns show separately the four dominant contributions to the 
systematic uncertainty: dependence on the B — > X c lv decay model, detector response, residual 
background, and B reco background subtraction. 



p* 

min 

(GeV) 


stat. sys. 


(M x 
X c Model 


) (GeV) 
Detector 


Background 


B reco Background 


0.9 


2.0728 ±0.0130 ±0.0127 


0.0021 


0.0092 


0.0077 


0.0037 


1.0 


2.0679 ±0.0120 ±0.0108 


0.0021 


0.0088 


0.0050 


0.0030 


1.1 


2.0609 ±0.0114 ±0.0097 


0.0023 


0.0087 


0.0033 


0.0016 


1.2 


2.0560 ±0.0109 ±0.0091 


0.0025 


0.0083 


0.0024 


0.0013 


1.3 


2.0465 ±0.0108 ±0.0086 


0.0024 


0.0080 


0.0015 


0.0013 


1.4 


2.0349 ±0.0112 ±0.0086 


0.0020 


0.0081 


0.0015 


0.0017 


1.5 


2.0344 ±0.0119 ±0.0089 


0.0022 


0.0081 


0.0014 


0.0026 


1.6 


2.0259 ±0.0133 ±0.0111 


0.0022 


0.0102 


0.0019 


0.0035 
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Table 2: Preliminary results for (M x ) for different values o/p^ in , with statistical and the total 
systematic errors. The last four columns show separately the four dominant contributions to the 
systematic uncertainty: dependence on the B — ► X c £u decay model, detector response, residual 
background, and B reco background subtraction. 



P min ( GeV ) 


stat. sys. 


(Mj 
X c Model 


k) ( GeV ) 
Detector 


Background 


B reco Background 


0.9 


4.366 ±0.049 ±0.057 


0.009 


0.034 


0.039 


0.023 


1.0 


4.338 ±0.043 ±0.046 


0.009 


0.033 


0.025 


0.016 


1.1 


4.300 ±0.040 ±0.038 


0.010 


0.032 


0.016 


0.006 


1.2 


4.276 ±0.037 ±0.035 


0.011 


0.030 


0.011 


0.006 


1.3 


4.239 ±0.036 ±0.031 


0.010 


0.028 


0.007 


0.006 


1.4 


4.190 ±0.036 ±0.029 


0.008 


0.027 


0.005 


0.007 


1.5 


4.180 ±0.038 ±0.028 


0.008 


0.026 


0.005 


0.006 


1.6 


4.146 ±0.042 ±0.034 


0.007 


0.031 


0.007 


0.009 



5 Systematic Uncertainties and Cross Checks 

Extensive studies have been performed to assess the systematic uncertainties in the measurement 
of (Mx) and {M x ) as a function of p^in- The ma i n sources of systematic errors are the subtraction 
of the combinatorical background of the B reco sample, the residual background estimate, inaccu- 
racies in the modeling of the detector efficiency and particle reconstruction, and the impact of the 
modeling of the hadronic states on the efficiencies and mass measurements. The uncertainty re- 
lated to the detector modeling and reconstruction is the dominant systematic error. Studies involve 
changes in the event selection and variations of the various corrections for particle reconstruction 
and resolutions, background suppression, and comparisons of results for various subsamples. The 
systematic errors are listed in Tables 1 and 2. 

The combinatoric B reco background is the principal background source. It is determined from 
the tues distribution as a function of Mx- The size of this background increases significantly for 
events with lower lepton momentum. The uncertainty in this background is estimated by varying 
the lower limit of the signal region in the rriES distribution. 

The residual backgrounds are due to misidentified hadrons and non-prompt leptons from 
decays and from semileptonic decays of and D s mesons produced via fragmentation. 
The processes leading to "right-sign" background from and D s decays contribute primarily 
to higher masses Mx- These decay processes dominate the background uncertainties since the 
combined uncertainty of their branching ratios is typically 30% [11]. Hadron misidentification 
contributes both "right-sign" and "wrong-sign" background. From a comparison of MC and data 
samples of events with "wrong-sign" leptons we conclude that the MC estimate of the hadron 
misidentification agrees with the data, and we take the 5% statistical error of this comparison as 
an estimate of the systematic error. 

The uncertainties in the detector performance and particle reconstruction have been estimated 
by Monte Carlo simulations of the track and photon efficiencies. Resolutions, fake rates, and 
background rates have been studied in detail by varying the adjustments to the MC simulation that 
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are introduced to improve the agreement with data. These include in particular the requirements 
on LAT and the minimum angular separation between a shower and the nearest track in the EMC. 
The root mean square (RMS) of the variation of the measured moments as a function of the lateral 
shower shape and the separation of the shower from the nearest charged track is assigned as the 
systematic error and is attributed to uncertainties in the simulation of hadronic interactions in 
the calorimeter. An error of 3% in the energy response of the EMC was assumed to estimate the 
uncertainty due to the photon energy measurement. The impact of track losses and fake tracks is 
largely independent of the minimum lepton momentum. 

The uncertainty in the acceptance variation is studied by varying all the branching ratios, in 
particular the composition of the hadronic high mass resonant and non-resonant states in the MC 
model. The observed bias is very small and we assign the RMS of the bias as the systematic error 
of this correction. The systematic error due to these variations is roughly constant for all p* n j„- 

The uncertainty in the Mx calibration is estimated by replacing the nine calibration curves by 
a single one. The observed change in the moments is assigned as the error in the Mx calibration. 
To verify the accuracy of the mass calibation method we apply it to a sample of B° — ► D* + £~u 
decays that are identified by partial reconstruction of the decay D* + — > D°ir + via the identification 
of the low momentum ir + . Figure 6 shows the hadronic mass-squared distribution for these events 
together with the measured {M x ) moments which agree well with the expectation to within the 
statistical errors and show a very minor dependence on p* min . 

The calibration method has also been extensively validated on MC samples. Since the moment 
measurements for different lepton momentum thresholds p^in are highly correlated, the moments 
were determined from the measured Mx and M x distributions for four separate intervals in the 
lepton momentum and compared with the true moments of MC generated distributions. The 
measured values agree well with the true values. 

A variety of tests are performed to test the stability of the moment measurements. The data 
are divided into several independent subsamples and the whole analysis is applied separately to 
each of them, including B reco background subtraction and mass calibration. The consistency of 
the measurement is checked for and B°, decays to electrons and muons, events from different 
run periods, with different E m i ss — \p m iss\, and of different B reco sample purity V . The results 
obtained from the two E m i ss — \p m i ss \ regions specifically address our understanding of the detector 
performance as this variable is very sensitive to the reconstruction of particles. A comparison of 
these results is shown in Figure 7. No significant systematic variations are observed. 
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Figure 6: Verification of the measurement of mass distribution and (M x ) moment for B — > D*£v. 
Left: The data (points) are compared to the MC prediction (histogram) for all selected events. The 
MC estimated background is shown as the shaded area at the bottom. Right: The measured {M x } 
moments ( data points ) as a function of p* nin are compared to the expectation, M^* , marked as a 
solid line. 
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Figure 7: A comparison of the results for (Mx) and (M x ) obtained from pairwise independent 
subsamples (statistical errors only). 
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6 Conclusions and Interpretation of Results 



We have performed a preliminary measurement of the first and second moment (Mx) and {M x } of 
the hadronic mass distribution in semileptonic B decays. These measurements have been carried 
out as a function of the threshold lepton momentum, p* min , extending from 0.9 GeV to 1.6 GeV. 
They are sensitive to the production of higher-mass charm states, both resonant and non-resonant. 
The reconstruction of a hadronic decay of one of the two B mesons and the kinematic fit to the 
full event have resulted in measurements with comparable statistical and systematic uncertainties. 
In addition, we have evaluated and taken into account the correlation between the individual data 
points. 

Figure 8a shows the moments (M x ) as a function oi p* min . As a result of several changes to 
the analysis and data selection, these new measurements fall below the preliminary measurements 
presented previously [1]. 

If this change were of purely statistical origin the level of compatibility of the two measurements 
would be approximately 5%. However, there is evidence that the systematic uncertainty related to 
the track and neutral selection used in [1] was underestimated. Since then, the total BB sample 
has doubled and the event selection, in particular the rejection of background tracks and photons, 
has been significantly improved. Furthermore, the moments are now extracted directly from the 
M x distribution, resulting in measurements that are much less sensitive to the uncertainties in the 
branching ratios and mass distributions of the individual hadronic states. 

The hadronic mass moments can be calculated in the framework of Heavy Quark Expansion 
(HQE), which expresses perturbative and non-perturbative corrections to the parton level calcula- 
tions in powers of the strong coupling constant, a s (rn,b), and in inverse powers of the B meson mass, 
1/rriB- We derive the two leading HQE parameters, Ai and A, that determine the non-perturbative 
corrections at 0{l/m 2 B ) from a minimum x 2 fit of the (M x ) moments measured for seven different 
lepton momentum thresholds. In these fits we take into account the strong correlations between 
the measured moments for different p* min - In these expansions we use the coefficients calculated 
in the MS mass scheme [3]. We use A 2 = 0.128 ± 0.010 GeV 2 as determined from the B - B* 
and D — D* mass splitting [12]. The six non-perturbative parameters pi 2 and 71,2,3,4 that appear 
at order C(l/m|), are expected to be of the order Aq CD and from dimensional arguments are 
assumed to be smaller than (0.5) 3 GeV 3 [12]. For this fit, we set their values to 71 = 0.0 GeV 3 , 
pi = ^(0.5) 3 GeV 3 , and we eliminate P2 by making use of the relation between p2,T~2 and Z4 and 
the D* — D and B* — B mass splittings represented in [12]. The fitted parameters result in a good 
description of the BABAR data, as indicated by the solid line in Figure 8. The fit result also agrees 
well with the measurement by the CLEO Collaboration [2] for lepton momenta above 1.5 GeV and 
is consistent with the DELPHI Collaboration [13] measurement corresponding to p^m = GeV 
(see Table 3). 

If we combine the first moment of the b — > S7 photon energy spectrum measured by CLEO [14], 
(J5 7 ) = 2.346 ± 0.034 GeV, with the CLEO (M%) measurement for p* min = 1.5 GeV and calculate 
Ai and A, the result does not give a consistent description of all available (M x ) measurements (see 
dashed line in Figure 8a). Here we used the coefficients for the HQE of (Ef) from [4]. The same 
effect is also illustrated in Figure 8b, where the result of the fit to the BABAR {M x ) measurements is 
displayed in form of a Ax 2 = 1 contour (39% confidence level) in the (Ai , A MS ) plane. The result- 
ing values for the HQE parameters are A MS = 0.53±0.09 GeV and Ai = -0.36±0.09 GeV 2 . These 
compare to earlier CLEO measurements of A MS = 0.35±0.07 GeV and Ai = -0.24±0.07 GeV 2 [2]. 
In Figure 8b, we also show the constraints from the hadronic mass moment measurements by CLEO 
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Table 3: Moment measurements from other experiments that are used in the fits. 



Experiment 


Quantity 


Value 




Hadron Mass Moments 


CLEO [2] 


{Ml) 


4.152 ± 0.066 GeV 2 


J 7 1 J 1 J I ILL 


{Ml) 


4.430 ± 0.080 GeV 2 




Lepton Energy Moments 


CLEO [15] 


R o 


0.6184 ± 0.0023 GeV 


CLEO [15] 


Rf 


1.7817 ±0.0013 GeV 


CLEO [15] 




0.6189 ± 0.0030 GeV 


CLEO [15] 




1.7802 ±0.0016 GeV 


DELPHI [13] 


Ri 


1.383 ±0.014 GeV 



Table 4: Measured quantities used in the determination ofTgL- 



Experiment 


Quantity Value 


BABAR [16] 
BABAR [17] 
BABAR [18] 
BABAR [18] 

HFAG 1 
BABAR [19] 
BABAR [20] 


t b o 1.523 ± 0.033 ps 
t b o 1.529 ± 0.031 ps 
t b o 1.546 ± 0.039 ps 
t b + 1.673 ± 0.039 ps 
f+/fo 1.053 ± 0.055 
B{B^Xiv) 0.1087 ±0.0035 
B(B -> X u lv) (2.14 ± 0.54) x 10~ 3 



and DELPHI, and from the {E^) moment from CLEO. The results of the {Ml) measurements are 
in good agreement. It should be emphasized that the one-o" areas indicated here do not include 
uncertainties due to the order 0(l/m B ) terms. 

Since the HQEs for the second hadronic mass moments are analogues to the total semileptonic 
rate, and since they contain the same non-perturbative parameters we can combine the seven 
measured moments with the semileptonic decay rate to perform a more comprehensive test of the 
HQE and to extract \V c b\ and the 6-quark mass, m&. For this extraction we choose the IS mass 
scheme as recommended in [4] because it leads to a better convergence of the perturbative expansion 
than the MS scheme. In this scheme, A 1S is related to the 6-quark mass via A 1S = m r /2-ml s . The 
coefficients for the expansions are taken from Bauer et al. [4]. We take into account the uncertainties 
of the other parameters of the expansion, i.e., for the independent 0(l/m B ) matrix elements we 
equate p\ = ^(0.5) 3 ± ^(0.5) 3 GeV 3 , % = 0.0 ± (0.5) 3 GeV 3 , and make the same assumption for p 2 
as stated above. For the semileptonic decay width, we take TgL = (4.37 ± 0.18) x 10 -11 MeV, as 
derived from individual BABAR measurements (see Table 4). 

Using the HQE expansions for the moments {Ml) and Tsl we determine three parameters by 
a x 2 minimization: \V c b\, m\ S , and Ax. The fit takes into account the experimental errors and 
their correlations. For treatment of the theoretical errors on the perturbative terms we follow the 

1 Heavy Flavor Averaging Group: /+//o from [11] has been rescaled by the lifetime ratio of B + to B° from [21]. 
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suggestions by Bauer et al. [4] to include errors estimated from dimensional analysis (dim) and 
from the size of the BLM term in the perturbative series [22]. These errors are added in quadrature 
and assumed to be 100% correlated for the various expansions. In addition, the sensitivity to the 
\jm\ terms is estimated from the change in the results of the x 2 minimization for a variation of 
the \jm\ matrix elements, assuming a flat, rather than Gaussian, error distribution. From these 
fits to the HQE expansions [4] for the BABAR hadronic moments and the semileptonic decay rate, 
based on BABAR measurements alone, we obtain 

\V cb \ = (42.10 ± 1.04(exp) ± 0.b2(dim © BLM) ± 0.50(l/?n|)) x 10~ 3 , 
m\ s = 4.638 ± 0.094(exp) ±0.062(dtmffi BLM) ± 0.065(l/m|) GeV, and 
Ai = -0.26 ± 0.06(exp) ± 0M(dim © BLM) ± 0.04(l/m|) GeV 2 . 

These results are in good agreement with other determinations [15, 19, 23]. The correlation between 
Ai and m\ s is 0.92, between \V c b\ and ml and Ai is —0.60 and 0.51, respectively. 

The shape of the lepton energy spectrum can be sampled in terms of moments that are also 
predicted by HQE calculations. Thus, to further test the HQE calculations we perform a separate 
fit using the measured lepton energy moments, as listed in Table 3, and the semileptonic rate T$l- 
Again, we are relying on the calculations of the expansion coefficients by Bauer et al. [4]. Since 
the x 2 minimization requires a proper treatment of the correlations among the measurements, we 
only consider published results that are either uncorrelated with all other measurements or for 
which the full covariance matrix is available. Thus we select four lepton moment measurements 
by the CLEO Collaboration [15] and one measurement by the DELPHI Collaboration [13]. In 
Figure 9(a,b) the fit results are shown separately for hadron and lepton energy moments, including 
the Ax 2 = 1 contours in the (\V c t,\ , m^ s ) and (Ai , m\ s ) planes, taking into account experimental 
and estimated theoretical errors (dim® BLM only) as explained above. The fit contours of hadronic 
and lepton energy moments do not overlap in either of the planes, but they each overlap with the 
one-cx constraint from the photon energy moments measured by CLEO. 

In conclusion, the measurements of the hadronic mass moments in semileptonic B decays pre- 
sented here agree well with measurements by other experiments of the same quantity, but have 
been extracted by a new technique which reduces the model dependence. By combining the second 
moments of the hadron mass distributions with earlier BABAR measurements of the semileptonic 
decay rate, we have been able to extract the HQE parameters m\ s and Ai and reduce the theoret- 
ical uncertainties on |V^,|. To further improve these measurements and examine the consistency of 
these results with measurements of other inclusive distributions in the context of HQE calculations, 
it will be very important to extend the currently available measurements to higher precision and 
thereby provide more stringent constraints on the 0(\/m? B ) parameters. 
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Figure 8: Comparison of results obtained from hadronic mass moments from the BABAR , CLEO 
and DELPHI experiments: a) the measured hadronic mass moments {M\) as a function of the 
lepton momentum threshold p* min - The solid curve is a fit to the BABAR data; the dashed curve is 
the HQE prediction obtained from measurements by CLEO of (E~) and (M^) for p* m in = GeV. 
b) Constraints on the HQE parameters A MS and \\ from: the fit to the BABAR (Mj^) measurements 
(Ax 2 = 1 ellipse), the CLEO and DELPHI (M^) measurements (one-a bands), and the CLEO 
b — > 57 energy moment (one-a band). The constraints are derived in the MS scheme to 0(l/m\). 
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Figure 9: Constraints on HQE parameters obtained from hadronic moments as measured by BABAR, 
all hadronic moments combined (BABAR, CLEO and DELPHI), and the combined lepton moments 
(CLEO and DELPHI), projected a) in the (\V c f,\, m\ s ) plane and b) the (\\, m\ ) plane. Also 
indicated is the constraint from the b — ► 57 photon energy moment measured by CLEO. 
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